Pulsar interpretation for the AMS-02 result 
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The AMS-02 collaboration has just published a high precision measurement of the cosmic positron 
fraction e^ l(e~ + e^), which rises with energy from ~ 5 GeV to ~ 350 GeV. The result indicates 
the existence of primary electron/positron sources to account for the positron excess. In this work, 
we investigate the possibility that the nearby mature pulsars are the primary positron sources. By 
fitting the data we find that the positrons from a single nearby pulsar, such as Geminga or Monogem, 
with the spectral index a ~ 2 can interpret the AMS-02 result. We also investigate the possibility 
that high energy positrons are generated by multiple known pulsars in the ATNF catalogue. Such 
a scenario can also fit the AMS-02 data well. Future precise measurements of fine structures in the 
positron spectrum would be a support to the pulsar scenario. 

PACS numbers: 96.50.S-, 97.60.Gb, 98.70.Sa 
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I. INTRODUCTION 

The Alpha Magnetic Spectrometer (AMS-02) collab- 
oration has just published its first result about the 
positron fraction in cosmic rays (CRs) with extremely 
high precision [l[. The positron fraction rises from ~ 5 
GeV continuously up to ^ 350 GeV, while the slope be- 
comes flat above ^ 100 GeV. This result is consistent 
with the PAMELA result about the positron fraction 
[2, |3| . Many studies show that primary electron/positron 
sources beyond the conventional cosmic ray (CR) model 
are necessary to explain the PAMELA data (see e.g. 
[J|). Astrophysical sources, like pulsars and pulsar wind 
nebulae (PWNe) |5l-[l7|. or dark matter (DM) annihi- 
lation/decay have been widely studied as the primary 
positron sources (e.g., |18j). 

Along the same line several works have appeared to ex- 
plain the AMS-02 data 0,[23. DM is still an attractive 
interpretation. As shown in Ref. [ly|, DM annihilates 
into T^T~ final states which results in a soft positron 
spectrum can account for the AMS-02 data quite well. 
Other interpretations that DM annihilation/decay into 
multiple /i or r leptons may also be fine to reproduce 
AMS-02 data. The DM annihilation scenarios require 
a large boost factor and the "Icptophilic" property of 
DM particles. However, significant secondary gamma 
rays are induced by cascade decay, final state radiation 
and inverse Compton scattering. Therefore, it is strongly 
constrained by the Fermi-LA T g amma ray observations 
from the Galactic center [2l|, IS^] or from dwarf galaxies 

[Hil. 

Pulsars are known to be powerful sources in the Galaxy 
to produce high energy electrons/positrons with energies 
of TeV scale or above |25l - [28| . Primary electrons are ex- 
tracted from the surface of the pulsar and are accelerated 
in the magnetosphere by strong electric fields (e.g. 10^^ 
V or higher) . Energetic curvature radiation is emitted in 
the strong magnetic field, which will result in electron- 
positron pairs due to the interactions with magnetic fields 
or low energy photons. The high energy gamma-ray pho- 
tons induced by those electron-positron pairs from pul- 



sars have been observed by Fcrmi-LAT |29j . 

Unlike the contributions from DM which are assumed 
to be continuously distributed in the halo and indepen- 
dent of time, the positron injections from pulsars are 
discrete in the Galactic disk and time-dependent. Since 
electrons/positrons lose energy quickly via synchrotron 
radiation and inverse Compton scattering when propa- 
gating in the Galaxy, the observed electrons/positrons 
above ^ 100 GeV can only come from a small range 
within a few kpc. A few nearby mature pulsars may 
have very important contributions to the high energy 
electron/positron spectrum and may induce significant 
deviation from the scenarios with continuous and steady 
injection. 

In this work, we investigate the possibility that nearby 
mature pulsars are the sources of the observed high en- 
ergy positions. We use a Markov Chain Monte Carlo 
(MCMC) method to fit the AMS-02 positron fraction 
data and determine the model parameters Q. To de- 
termine the properties of the electron backgrounds, we 
also include the PAMELA electron data |30l in the fit. 
We consider the possibilities that a single nearby pulsar 
such as Geminga or Monogem is the source to produce 
the observed positrons. Through fitting to the data we 
get the constraints on the parameters of a single pulsar, 
such as the distance, age and total injected e^ energy. 
High energy e may also be generated by multiple pul- 
sars rather than a certain single pulsar. We then discuss 
the positron spectrum from a population of pulsars based 
on the ATNF pulsar catalogue [31| . The multiple pulsars 
may produce bump like structures in the positron spec- 
trum. We discuss the possibility to distinguish the pul- 
sar and the DM scenarios by a future experiment, such 
as the Chinese satellite experiment DArk Matter Particle 
Explorer (DAMPE), which is planned to be launched in 

2015 m. 

It is worth emphasizing that the number and the en- 
ergy distribution of e pairs injected from the pulsar 
magnetosphere are still open questions. Since there may 
exist a PWN between the pulsar and the interstellar 
medium (ISM) (or a supernova remnant), the spectrum 



of e pairs from the pulsar magnctosphcre would be mod- 
ified by the termination shock and radiation cooling be- 
fore they are injected into the ISM [la. Il7l|. In compar- 
ison with the DM scenario, it is very difficult to obtain 
a concrete form of the initial e spectrum from pulsar 
models. Therefore, there would be large uncertainties in 
the pulsar scenario to explain the high energy positron 
excesses. 

The paper is organized as follows. In Sec. II we de- 
scribe our treatments for the CR backgrounds and prop- 
agation parameters. In Sec. Ill, we discuss the injec- 
tion e^ spectra from the pulsar and their propagations. 
In Sec. IV, we consider the possibility of a single pul- 
sar as the source of the observed high energy positrons, 
and take Gcminga and Monogem as benchmark exam- 
ples. In Sec. V, we calculate the positron contributions 
from multiple pulsars. Then we discuss the possibility 
to distinguish the pulsar scenario from the DM scenario 
by a future experiment in Sec VI. Finally Sec. VII is our 
conclusion and discussion. 



II. COSMIC RAY BACKGROUNDS 



account the solar modulation effect [40[ • Note the low en- 
ergy part (< 5 GcV) of the positron fraction may not be 
easily explained with the simple solar modulation model, 
and more complicated charge-sign dependent modulation 
is necessary [4l|, \^ . 



III. HIGH ENERGY e± PAIRS FROM THE 
PULSAR 

The energy of e* injected into the ISM is limited by 
the rotational energy loss rate of the pulsar. The rota- 
tional frequency il = 2Tr/P, with P the pulse period, de- 
creases as f2 = —ail". Here n = flfl/Cl^ is the breaking 
index which can be calculated from the measurements 
of 57, 51 and tl. The upper limit of the total e^ en- 
ergy is determined by the pulsar spin-down luminosity 
E = m\fl\ = af7"+\ where / = (2/5)M,i?2 js the mo- 
ment of inertia of the pulsar, Ms and Rg are the mass 
and radius of the pulsar, respectively. For the magnetic 
dipole radiation the braking index n = 3, and the ro- 
tational frequency fl is given by (see Ref. 0, Q and 
references therein) 



The background to explain data includes primary 
electrons from the CR sources, and the secondary 
positrons/electrons which are generated in the collisions 
between the CR nuclei and the ISM. In this work, we use 
the GALPROP code to calculate the background JHHI. 
We employ the diffusion reacceleration model for CR 
propagation. The propagation parameters are adopted 
by fitting to the Boron-to-Carbon ratio and unstable-to- 
stable Beryllium ratio [3a] (see also (Sg]). The parame- 

5 = 0.377, 



ters are £'oUo=4GV 

va = 36.4 km s^^ and Zh 



5.94 X 1028 cm^ s 



4.04 kpc. With these propa- 
gation parameters, the injection spectrum of the protons 
is fitted according to the PAMELA [s^l and CREAM 
[3a | data. The fitting parameters of the proton injection 
spectrum are vi = 1.80, 1^2 = 2.36 and R^^ — 11.7 GV, 
where vi and V2 are the spectral indices below and above 
the break rigidity R^^ ]1^ . 

The secondary positrons and electrons can be calcu- 
lated according to the pro ton spectrum and the propa- 
gation model (see e.g. |39[). To involve some unknown 
uncertainties, e.g., from the ISM density distribution, the 
hadronic interactions and the nuclear enhancement factor 
from the heavy elements, we rcscalc the calculated fluxes 
of secondary electrons and positrons with a free factor 
Cg+ in order to fit the data. The injection spectrum 
of the primary electrons is parameterized by a broken 
power-law with respect to the rigidity (or momentum), 
q{R) oc {R/R^^)~''^^''^, with 71 and 72 the indices below 



and above the break rigidity R^^. A further normaliza- 
tion factor Ae is needed in the fitting procedure. 

For energies less than several tens of GeV, the fluxes 
of CR particles will be modulated by the solar environ- 
ment, known as solar modulation. The force fleld approx- 
imation, with only one parameter 0, is used to take into 
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where f2o is the initial rotational frequency of the pulsar. 
To = Sc^I/B^R^riQ [Bg is the surface magnetic field) 
is a characteristic time scale describing the spin-down 
luminosity decays, tq is usually assumed to be ^ 10^ 
year [27| . The spin-down luminosity of a pulsar is then 
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For the pulsar with t :^ tq, E decreases as t^^. The age 
of the pulsar T can be obtained through the rotational 
energy loss rate approximately, T = 51/2 |f2|. The total 
energy of electrons and positrons injected from a pulsar 
is assumed to be proportional to the rotational energy 
loss, which is 



• T2 
Eout = »7e± / Edt ~ r]a±E — , 

To 



(3) 



where ri^± is the fraction of the rotational energy con- 
verted into the energy of electrons and positrons. 

The propagation of high energy e can be described 
by the diffusion equation in the spherically symmetric 
approximatio: 
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where f(r,t,E) is the time dependent differential den- 
sity of electrons and positrons (the differential fiux is 
cf{r,t,E)/47r), D{E) oc I3Do{E/EqY is the diffusion 
coefficient with Dq and 5 the same as the background 



calculation (Sec. II). The energy loss rate due to syn- 
chrotron and inverse Conipton scattering is adopted as 
b{E) = boE"^ with 60 = 1-4 X 10-16 GcV-i s"! [ll|. 

For the burst-like source, the source term can be taken 
as a (5 function Q{r,t,E) ex 5{r — ro)S{t — to). This is 
a good approximation for the pulsar with T ^ tq since 
most of the rotational energy is lost during the time scale 
tq. The injection energy spectrum of the pulsar is pa- 
rameterized as a power-law function with an exponential 
cutoff 

QiE,r,t) ^ QoE-'^ cxp{~E/E,^,)S{r~ro)S{t^to), (5) 

where a is the spectral index, Ecut is the cutoff energy, 
Qo is the normalization factor which can be determined 
by the total injection energy E^ut- The solution of Eq. 
dlD for source term Eq. ([5|) is f26l. [27} 
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where d is the distance of the pulsar from the earth and 
td is the diffusion time into the ISM. Note that td may be 
different from the actual age of the pulsar T since elec- 
trons and positrons may spend some time in the PWN 
before their injection in the ISM. Here we simply assume 
td « T. E'max — ibotd)~^ is the maximum energy of elec- 
trons and positrons surviving from cooling. For e^ with 
energies larger than E'max, f{d,td,E) is taken to be 0. 
The diffusion distance r^n is given by 



rdit{td,E)^2jD{E)td 
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Note that for old pulsars with T ^ 10'' yr, the ( 
tion energy rate is suppressed by l/T^. In addition, the 
positrons from old pulsars should not contribute much to 
the observed flux at high energies due to the energy loss 
in the ISM. The upper limit of the propagation time of 
the e* with certain energy is 



1 
b^ 



t < T^ - 2.3 X 10^ yr 



VTevj 



(8) 



For young pulsars with T < O(IO^) yr, the positrons 
may not have enough time to propagate to the Earth. 
Moreover, these positrons may still be trapped in the 
PWN and not injected in the ISM during such short time 
scale. The limit of the propagation time of the e* also 
suggests an upper limit of the diffusion distance [ifl] , 



<2^/D{E)t^ Ikpc 
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(9) 



Thus, the nearby pulsars with ages T ^ O(IO^) yr and 
distances d ^ 1 kpc are thought to be good candidates 
to interpret the exotic high energy positrons. 



IV. SINGLE PULSAR 

In this section, we study the possibility to use a single 
pulsar as the high energy positron source to fit the AMS- 
02 positron fraction. The background electron spectrum 
is described by four parameters, i.e. Ae, 71, 72, R^^, 
Ce+ and 4>e which are free parameters in our fit. The 
positron/electron spectrum from pulsars are described 
by five parameters in Eq. ©, the distance d, the prop- 
agation time of <d, the total injected energy i?out, the 
cutoff energy i^cut and the spectral index a. For nearby 
known pulsars, the distance and age are adopted by the 
ATNF catalogue data. Note that the measurements of 
the pulsar distance still have some uncertainties. More- 
over, the fact that pulsars have velocities of ~ O(IO^) 
km s"^ means that the current distance of the pulsar is 
different from that during the e injection period (see 
e.g. Q and references therein). Here we do not take into 
account such uncertainties in our fit. 

We consider two nearby pulsars Geminga 
[J0633+1746] with a distance of d = 0.25 kpc and 
age of T = 3.7 x 10^ yr, and Monogem [B0656+14] 
with d = 0.28 kpc and T = 1.1 x 10^ yr. Here we take 
the age of the pulsar as the diffusion time. We use a 
MCMC method to determine -Ecut, a and i^tot- The 
resulting positron fraction and electron spectrum for the 
best fitting parameters for Geminga and Monogem are 
shown in Fig. [TJ The best fitting parameters are given in 
Table m From Fig. [T]we can see that the pulsar scenario 
with these parameters can fit the AMS-02 positron 
fraction and the PAMELA electron flux data very well^ . 
Since Monogem is younger than Geminga, high energy 
positrons from Monogem would have smaller energy 
loss. Thus, the cutoff of the Monogem spectrum is larger 
than that of Geminga. It is possible to observe such 
cutoff with more AMS-02 data accumulation. 



pulsar 



Geminga 
Monogem 



E, 



a Eu 



1.0 
0.62 



1.98 
2.04 



14.2 
3.30 



log(yle) 71 72 Rl 



Ce+ 



-8.93 1.74 2.75 3.61 1.53 720 
-8.93 1.75 2.75 3.62 1.61 735 



TABLE I: Parameters of best fit for two nearby pulsars (Scut , 
a and Etot) Geminga (with d = 0.25 kpc and T = 3.7 x 10^ 
yr), Monogem (with d = 0.28 kpc and T = 1.1 x 10^ yr), and 
electron backgrounds (Ae, 71, 72, Rl,^, c^+ and 0e). -Rbn -E-cut 
and -Bout are in units of MV, TeV and 10*^ erg, respectively. 



Aa is normalized at 1 MeV in unit of cm 



MeV" 



We also show the confidence regions in the plane of 



^ Note that here we have considered the low energy electron data 
from 1 GeV to 20 GeV in our fit. Since this setting tends to 
choose a softer spectrum for the electron background, our fit 
shows a tension between PAMELA and Fermi/HESS electron- 
positron data. More discussions about this issue can bo found in 
Rofs. [li,[5i. 
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FIG. 1: The positron fraction (left) and electron flux (right) for the exotic e from Geminga [J0633+1746] (upper) and 
Monogem [B0656+14] (lower). Also shown are the positron fraction data from AMSOl 0iHEAT94+95 Q, HEATOO ^, 
PAMELA 01 and AMS02 0], and electron flux data from PAMELA Q, ATIC [H], HESS [13, [3 and Fermi-LAT ^. 
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FIG. 2: la and 2a confidence regions in the plane of a vs. 
Etot for Geminga (blue) and Monogem (red), respectively. 
Horizontal lines denote the maximum Eout derived by Eq. 



Q with TO 
best fits. 



10 yr and ri^± — 1. Also shown are points for 



a vs. Etot in Fig. [2] Compared with the parameters 
required to fit the previous PAMELA positron fraction 
data, the spectrum of the single pulsar becomes softer. 
The typical spectral indices for Geminga to interpret 
the AMS-02 data is 1.8 -- 2.1, while it is -- 1.5 to fit 
PAMELA positron fraction data given by Ref. [5|. The 
observed spin-down luminosity of Geminga is 3.2 x 10'^^ 



erg s~"^. By using Eq. ([3]) and adopting tq = 10^ yr, the 
total injection energy to electron and positron pairs is 
1.2?7g± X 10"*^ erg which is comparable with the typical 
values of our best fit. 

A similar conclusion can be applied for Monogem. The 
typical injection spectral of Monogem to interpret the 
AMS-02 data has a power-law index within 1.9 ^ 2.2 
which is softer than that for the PAMELA positron 
data. The observed spin-down luminosity of Monogem 
is 3.8 X 10'^^ erg s~^, resulting in the total injection en- 
ergy 1.4877e± X 10^* erg. The typical value of the E'tot in 
our fit is within 2 '--^ 6 x 10"*^ erg. 

To relax such energy tension, one can change propa- 
gation model. Another possibility is that Geminga or 
Monogem may not be the only source to contribute to 
the observed high energy positrons. 

Then we investigate the correlations between the pul- 
sar distance and age with fixed spectral index and the 
cutoff energy following [lg| . Applying an MCMC method 
to fit the positron fraction and electron spectrum with 
a = 2 and E^ut = 1 TeV, we derive the la and 2(t con- 
fidence regions in the pulsar parameter space, as shown 
in Fig. [3 

In Fig. El 177 selected pulsars in the ATNF catalogue 
with d < 3kpc and 5 x 10^ < T(yr) < lO'^ are also 
plotted, ii^tot for each pulsar is estimated from E and 
T by Eq. ([3]) with tq ~ 10"* yr and 77e± varying from 5% 
to 50%. From Fig. [3] we note that the 7 nearby pulsars 
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FIG. 3: la (solid contour) and 2a (dashed contour) confidence regions in the pulsar parameter space to fit the positron fi-action 
and electron spectrum for a — 2 and -Bcut ~ 1 TeV. The regions are projected onto the Etot vs. d (upper left), Etot vs. T 
(upper right), and T vs. d (lower) planes. The circles and bars denote the 177 pulsars in the ATNF catalogue within d < 3kpc 
and 5 x 10'' < r(yr) < 10^. Etot for each pulsar is estimated in the range of 5% < rj^± < 50%. 7 nearby pulsars with d < 0.5 kpc 
and T < 10® yr are marked by colors. 



with d < 0.5 kpc and T < 10^ yr are more likely to fit the 
data. We have marked them by colors, with the color 
green/magenta corresponding to Geminga/Monogem. 

Fig.[3]shows that the favored region in the pulsar space 
is rather small. There are several pulsars located near the 
favored region, especially 7 nearby pulsars marked by col- 
ors. These pulsars could also have sizable contributions 
to the high energy e flux. Therefore, a more reasonable 
treatment may include the contributions of all suitable 
pulsars. 



V. MULTIPLE PULSARS 

It is possible that the flux of high energy elec- 
tron/positrons arc contributed by many pulsars. There- 
fore we sum the contribution from all the 177 mature 
pulsars in the ATNF catalogue to get the positron flux 
following the method in Ref. [ll|. For each pulsar, we 
randomly assign the parameters in the following ranges: 
700 < £;cut(GcV) < 3000, 1.5 < a < 2.3 and 5% < 
Ve± !i 30%. The results arc shown in Fig. 2] Obviously, 
by summing the contributions of all pulsar, even low e 



pair conversion efficiency r]^± could be enough to fit the 
data. 

Since the pulsar parameters vary in large ranges, the 
resulting total spectrum of all pulsars also varies in a 
wide band. The energy cutoff of each pulsar depends on 
the minimum of the injection cutoff Ecut and the cooling 
cutoff -Emax, and is different from each other. It is further 
shown that few nearby pulsars could dominate the total 
flux. Therefore, for each combination of the parameters, 
the sum spectrum tends to have several bumps at high 
energies, as shown in Fig. [J] 

Since some pulsars radio beams do not point toward 
the earth, the ATNF catalogue is incomplete. There 
might be a diffuse population of pulsars which are be- 
yond the observed catalogue. This diffuse component 
may contribute as another "background" of the elec- 
trons/positrons. Similar as done in [l9|, we introduce a 
continuously distributed source component of the diffuse 
pulsars, with spatial distribution [51| 

Q{R, ^)°^ i^j exp 

(10) 
where Rq = 8.5 kpc and z^ ~ 0.2 kpc. The energy spec- 



r 5.56{R-Rq)] 


/ 


kl 


— 


cxp 




L ^0 J 


v 


^s 



10" ^ 



AMS ^^ PAMELA 
HEAT94+95 ' — ■ — ' AMS02 
HEATOO Bkg 



10'' 




ATIC 

HESS 

Fermi-LAT 



PAMELA 
Bkge^ 
Bkg e + e' 



E (GeV) 



E (GeV) 



FIG. 4: The contributions to the positron fraction (left) and electron flux (right) of all the 177 pulsars in the ATNF catalogue 
with d < 3kpc and 5 x 10^ < T(yr) < lO''. The parameters for each pulsar are randomly assigned in the following ranges: 
700 < -Ecut(GeV) < 3000, 1.5 < a < 2.3 and 5% < ri^± < 30%. For a particular combination of the parameters, the resulting 
spectrum including the contributions of all pulsars is represented by a grey line, while a golden line shows only the exotic 
contributions from the pulsars. A representative choice is shown by black lines (the solid line for the total result, and the 
dot-dashed line for the exotic contribution). 
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FIG. 5: Illustration of the positron fraction from both the dif- 
fuse pulsars and the 7 nearby mature pulsars. The parameters 
of the 7 pulsars are chosen randomly. 



truin of the diffuse pulsars is parameterized by a power- 
law function with an exponential cutoff. 

Fig. [5] shows an illustration of the CR positron frac- 
tion from both the diffuse pulsars and the 7 nearby ma- 
ture pulsars mentioned above. Here the energy spectrum 
of the diffuse pulsars is adopted to be proportional to 
E'"^ exp{—E /600 GeV). The parameters E'cut, a and r]e± 
of the 7 nearby pulsars are chosen randomly in the ranges 
described above. Note that only one combination of the 
parameters are shown in in Fig. [5] as an illustration. 



VI. DISTINGUISH PULSAR FROM THE DM 
SCENARIO 

Both the pulsar and the DM scenarios can fit the AMS- 
02 data well [19j . It is a fundamental problem to distin- 



guish these two scenarios. As we discussed above if the 
positron excess is from a few nearby pulsars, it may have 
a characteristic spectrum with many structures. If such 
fine structures are discovered, it would be a strong sup- 
port to the origin of multiple pulsars for the high energy 
electrons/positrons. Here we explore the possibility of 
distinguishing such two scenarios by using future elec- 
tron/positron spectrum observations. The similar dis- 
cussions can be found in Refs. [3, l9l.[la|. 

We generate the mock data using the pulsar set de- 
noted by black dot-dashed line in Fig. |4] as an exam- 
ple, and consider the observation capability of DAMPE. 
The mock data are produced following the method in 
Ref . [lal . The number of particles detected in a certain 
energy bin with a gaussian energy resolution is given by 



N{E) ^ AtSE A / dE'(j){E') 



V2^ 



: CXp 



(E' - Ef 



2a' 



(11) 

where (^ is the differential flux, 8E is the width of the 
energy bin, Ai is the observation time, A is the ge- 
ometrical factor of the detector and a is determined 
by CT = /S.E/2. AE is adopted to be the form of 
AE/E = a/y/E/GeY S) b, which is normalized to 1.5% 
and 10% at 1 TeV and 1 GeV for DAMPE, respectively 
[33j . The geometrical factor and performing time of the 
detector are taken to be 0.5 m^ sr and 5 yr, respec- 
tively [33. The relative statistical uncertainty can be 
estimated by ^/Ne/Ne- The systematic uncertainty is 
assumed to be mainly determined by the capability of 
distinguishing electrons/positrons and other CR parti- 
cles [131. Here the e/p separation is taken to be 5 x 10^ 
corresponding to the detector thickness of 32 radiation 
lengths [33] ■ The relative systematic uncertainty is esti- 
mated by {Np/5 X 10^)/Ne. 

We use the MCMC method to explore the possible 
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FIG. 6: The positron fraction (left) and electron flux (right) for the background together with an exotic component from 
multiple pulsars or DM annihilation in t'^t^ channel. The "mock" DAMPE data are assumed to be contributed by multiple 
pulsars and are generated using the pulsar set denoted by black dot-dashed line in Fig. [l] The DM contribution corresponds 
to the best fit of an MCMC parameter scan. 



DM parameter space to fit the mock electron/positron 
flux data for DAMPE, the positron fraction data from 
AMS-02 and the electron flux data from PAMELA. As 
above, the proton injection spectrum is fixed, while the 
parameters of the primary electron injection spectrum 
arc free. Therefore, the background corresponding to the 
fit for the DM is usually different with that for multiple 
pulsars. The DM annihilation final states are taken to 
be T~^T~ , which induce a soft positron spectrum favored 
by the AMS-02 result. From the results shown in Fig. IH 
we find that the behavior of the electron/postron spec- 
trum from the DM source is mainly determined by the 
mock DAMPE data below ~ 300 GeV due to very small 
uncertainties, and it cannot reproduce the fine structures 
above 300 GeV at the data. If the differences between 
the electron/positron spectra from the DM and the pul- 
sar origins are significant as the examples shown here, 
DAMPE would have the capability to discriminate these 
two scenarios. 



sar derived from Eq. ([3]). Considering the uncertainties 
from CR propagation parameters and the pulsar models, 
such a tension may be relaxed. We then consider the 
case that the positrons are from multiple pulsars in the 
ATNF catalogue. We find such scenario can also fit the 
AMS-02 data very weU. 

It is shown that pulsars can be a natural explanation 
of the AMS-02 positron excess. Compared with the DM 
scenario, pulsar scenario may have some distinct features 
to be distinguished from the DM models. It is very 
possible that there might be fine structures on the elec- 
tron/positron spectrum in the pulsar scenario, because 
the parameters of pulsars might differ from one to an- 
other [9|. Furthermore, since one or several nearby pul- 
sar(s) may dominate the flux of high energy positrons, 
it may induce a remarkable anisotropy [5|, llq . Both the 
fine structures and the anisotropy could be tested with 
future observations. 
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